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ABSTRACT 

Four quiescent elliptical galaxies were imaged with the NRAO VLA at 8.5 GHz. 
\ Within the context of canonical advection-dominated accretion flows (ADAFs), these 

VLA images plus published black hole masses constrain the accretion rates to be < 

O ■ 1.6 X 10"^, < 3.6 X 10"'', < 7.8 x 10"^, and < 7.4 x 10"'' of the Eddington rates. These 

■ , _ , _ 

■ ADAF accretion rates derived at 8.5 GHz have important implications for the levels of 

■ soft and hard X-rays expected from these quiescent galaxies. 

O 

Q^i Subject headings: accretion, accretion disks - galaxies: individual (NGC 4291, NGC 4564, 

Q , NGC 4621, NGC 4660) - galaxies: nuclei - radio continuum: galaxies - X-rays: galaxies 



MOTIVATION 



Evidence is accumulating that nearby galactic nuclei commonly harbour massive dark objects 
5^ \ (Magorrian et al. 1998). These objects are probably black holes, because such remnants from the 

QSO era should be common (Fabian & Rees 1995) and because some nuclear star clusters with the 
requisite mass and size would be improbably younger that their host galaxies (Maoz 1998). Given 
a black hole mass, Mbh, a next important step is to constrain the rate, M, at which material 
is being accreted onto the black hole. Associated with any black hole is its Eddington accretion 
rate. Me, which is the rate necessary to sustain the Eddington luminosity. For a 10% radiative 
efficiency and Mbh = lO^mg M©, Me = 2.2m8 Mq yr~^. For giant elliptical galaxies hosting 
black holes, accretion-rate estimates have long been available for comparison with Eddington rates 
(Fabian & Canizares 1988; Mahadevan 1997). These estimates are based on Bondi (1952) accretion 
from an interstellar medium with temperature T = lO^Ty K and pressure P = 10^ Pg cm~'^ K, 
with the accretion being characterized by a Bondi radius r-Q ^ A.SmsTj pc and a Bondi rate 



^Current address; Renaissance Technologies Corporation, 600 Route 25A, East Setauket, NY 11733-1249 



- 2 - 



Mb ~ 1.9 X lO^^mgPeTV Ms yr~^- Then a black hole with = 1 in a medium with Pq = 1 
and T7 = 1 will accrete at a Bondi rate Mb ~ 1-9 x 10~^ Mq yr~^, which is almost four orders 
of magnitude less than the associated Eddington rate of Me = 2.2 Mq yr~^. The Bondi rate 
estimates are therefore extremely sub-Eddington. Still, such low rates could be pervasive among 
nearby ellipticals and, moreover, could define a minimum level of activity for ellipticals hosting 
massive black holes 

This Letter examines the consequences, in the radio and X-ray regimes, of such low accretion 
rates in nearby elliptical galaxies. The approach is to obtain deep radio continuum images of 
four ellipticals studied by Magorrian et al. (1998) and then interpret those images within the 
context of canonical advection-dominated accretion flows (ADAFs) , reviewed recently by Narayan, 
Mahadevan, & Quataert (1998). The radio continuum from an ADAF is thermal synchrotron 
emission from a magnetized plasma, and the conversion between the ADAF radio power. Pi,, and the 
ADAF accretion rate, Ma, is relatively straightforward if the black hole mass is known (Mahadevan 
1997; Yi &; Boughn 1998). Recent efforts along these lines have focused on galaxies with significant 
radio emission from nonthermal synchrotron jets (Di Matteo et al. 2000). In constrast, this study 
examines four quiescent elliptical galaxies previously undetected at radio wavelengths (Wrobel 
1991; Wrobel &: Heeschen 1991), thereby minimizing one of the largest potential complications in 
an ADAF analysis - the role of jet outflows. Further evidence for the quiesence of these four galaxies 
comes from their use as absorption-line templates in optical spectroscopic studies (Ho, Filippenko, 
Sz Sargent 1997), plus their weak or undetected X-ray emission (Beuing et al. 1999). For these four 
quiescent ellipticals, the new radio imaging constrains the ADAF accretion rates to be < 1.6 x 10~^, 
< 3.6 X 10""^, < 7.8 X 10""^, and < 7.4 x 10~^ of the Eddington rates. These ADAF accretion rates 
derived at 8.5 GHz have important implications for the levels of X-ray emission expected from these 
quiescent galaxies. 

2. OBSERVATIONS AND IMAGING 

The Very Large Array (VLA), described by Thompson et al. (1980), was used 1999 May 6-7 
UT in its 1-km configuration to observe NGC4291, NGC4564, NGC4621, and NGC4660. Data 
were acquired in dual circular polarizations and at a center frequency 8.4601 GHz with bandwidth 
100 MHz. Observations were made assuming a coordinate equinox of 2000. The phase calibrator 
J1153+8058 was used for NGC4291, while J1239+0730 was used for the other three galaxies. 
Scheduled on-source integration times were 105-120 minutes. The galaxies were observed at a 
priori positions from the Digitized Sky Survey (DSS) web site. Improved DSS positions with 1-D 
error estimates, (TdsSj have since become available (Cotton, Condon, &; Arbizzani 1999) and will be 
used below, cross depends in part on the galaxy major axis, ^m, and values for both terms appear 
in Table 1. Observations of 3C 286 were used to set the amplitude scale to an accuracy of about 
3%. The data were calibrated using the 1999 October 15 release of the NRAO AIPS software. No 
self-calibrations were performed. 
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The AIPS task IMAGR was used to form and deconvolve images of the Stokes IQU emission 
from each galaxy. All images were made with natural weighting to optimize sensitivity. The Stokes 
/ images appear in Figure 1, with crosses indicating the the DSS positions and their errors. The 
large ellipses shown in three panels outline the galaxies' major and minor diameters and elongation 
position angles (Wrobel 1991). The corresponding ellipse for NGC4621 falls beyond the panel 
boundaries. For each panel the table gives the image resolution, i?8.5GHz, expressed as the FWHM 
dimensions and elongation orientation of the elliptical-Gaussian restoring beam; the image rms 
value, (T8.5GHz, in units of microjanskys per beam area (^Jy ba^^); and the image scale, s, based on 
the distance, D, from Magorrian et al. (1998) for Hq = 80 km s"-*^ Mpc^^. Each galaxy image was 
searched for an emission peak >S'8.5GHz > 4a"8.5GHz within a circular region of radius 4 ctdss- The 
resulting peak strengths appear in the table. For the two detected galaxies, quadratic fits in the 
image plane yielded the peak flux densities tabulated, along with their errors that are quadratic 
sums of a 3% scale error and 0"8.5gHz- Those fits also yielded the following radio positions, along 
with their 1-D errors, itvlAj which are dominated by the peak signal-to-noise ratio (Ball 1975): for 
NGC4621, a(J2000) = 12'^42™02^50, (5(J2000) = 11°38'48".6, and (Jvla = 0".4; for NGC4660, 
a(J2000) = 12^* 44^32^35, (5(J2000) = ll°ll'26".6, and avLA = 0".7. Each detected galaxy is less 
than 40% linearly polarized. 



3. IMPLICATIONS 

Radio continuum can be due to thermal synchrotron emission from an ADAF onto a black hole, 
nonthermal synchrotron jets emerging from an ADAF, or nonthermal synchrotron jets emerging 
from a standard accretion disk which itself emits no radio photons. In the radio regime, this Letter 
analyzes only thermal synchrotron emission from an ADAF. The distances, D, from Magorrian 
et al. (1998) are used to convert the new 8.5-GHz detections or limits to the observed powers, 
-Ps.SGHz) given in the table. These rank among the the most sensitive radio powers, or power limits, 
available in the literature for elliptical galaxies (Wrobel & Heeschen 1991). Because emission from 
nonthermal synchrotron jets could account for some unknown fraction of the radio detections, the 
observed peak flux densities will in general impose upper limits to any ADAF powers, P8.5GHz- 

The radio power, P^, predicted for an ADAF at a frequency, depends upon: a, the standard 
Shakura-Sunyaev parametrization of the kinematic viscosity in the accretion flow; /?, the ratio of 
gas pressure to total pressure; ms, the black hole mass in units of 10^ Mq] rriA, the ADAF accre- 
tion rate normalized to the Eddington rate; and Tg, the equilibrium temperature of the electrons 
(Mahadevan 1997; Yi & Boughn 1998). Canonical ADAFs adopt a = 0.3 and an equipartition 
P = 0.5. Moreover, for most cases of interest Tg is fairly insensitive to the other model parameters, 
so an intermediate value of 2.5 x 10^ K is adopted. These values for a, and Tg then predict that, 
at frequency = 8.5 GHz, 

P8.5GHZ = 6.7 X lO^imfmi W Hz'^ (1) 
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The values for Ps.sGHz from this work and for ms from Magorrian et al. (1998), when inserted 
into Equation (1), result in the normalized ADAF accretion rates, rriA, given in the table. These 
normalized ADAF rates are highly sub-Eddington, being < 1.6 x 10""^, < 3.6 x 10"'', < 7.8 x 10""^, 
and < 7.4 x 10~^ of the Eddington rates. The latter two cases correspond to the galaxies with radio 
continuum detections. Canonical ADAFs exisit only at normalized accretion rates below a critical 
rate, merit ^ 0.3q;^ = 0.027 (Esin, McClintock, &; Narayan 1997). The normalized rates, niA, in the 
table easily meet this consistency check for canonical, and other reasonable, values of a. For the 
four quiescent ellipticals in this study, the black hole masses from Magorrian et al. (1998) yield the 
tabulated Eddington accretion rates, Mg. Those rates are used to convert the normalized ADAF 
rates, m\, in the table into the absolute ADAF rates. Ma, also entered in the table. 

The ADAF accretion rates derived at 8.5 GHz, in both normalized (ttia) and absolute (Ma) 
form, have consequences for the predicted levels of X-ray emission from these quiescent ellipti- 
cal galaxies, and these predicted levels must not violate the observed levels. Only one galaxy, 
NGC4291, is a weak X-ray emitter and the other three remain undetected in the ROSAT All Sky 
Survey (Beuing et al. 1999). That survey was conducted with the Position Sensitive Proportional 
Counter at soft X-rays (0.5-2 keV) and at an angular resolution of ~ 0.5'. The table gives values 
for, or upper limits to, the observed soft X-ray luminosities, LrasSi scaled to the Magorrian et al. 
(1998) distances. No deprojection analysis has yet been attempted for NGC4291, since no ROSAT 
data from the High Resolution Imager are available. Data at hard X-rays (2-10 keV) are also 
lacking. 

An elliptical galaxy harboring an ADAF will be a source of X-rays from the ADAF itself 
and from the galaxy's interstellar medium. Each photon source will be discussed in turn. The 
normalized ADAF accretion rates derived at 8.5 GHz are so low that bremsstrahlung emission, not 

inverse Compton scattering, will dominate the ADAF's X-rays (Mahadcvan 1997; Yi & Boughn 
1998). The latter authors provide an expression for the bremsstrahlung luminosity of a canonical 
ADAF that, at the Tg adopted for Equation (1), reduces to 

^A,ikeV = 1-2 X lO^^msmi erg s"^ (2) 

at a fiducial soft X-ray frequency of = 2.42 x 10^^ Hz; and reduces to 

^A,6keV = 7.2 x lO^Vgrhi erg s'^ (3) 

at a fiducial hard X-ray frequency of u = 1.45 x 10^^ Hz. Applying the tabulated values for mg 
and rriA to Equations (2) and (3) yields the tabulated predictions for the soft and hard X-ray 
luminosities of the ADAFs. The predicted luminosities at 1 keV are consistent with the published 
ROSAT hmits for NGC4564, NGC4621, and NGC4660, as well as for the ROSAT detection of 
NGC 4291 but only if that detection is dominated by the galaxy's interstellar medium rather than 
by its ADAF (Beuing et al. 1999). Note further that the ADAFs are expected to be six times more 
luminous at hard than at soft X-rays, so these galaxies must must be considered prime targets for 
observations in the 2-10 keV region with the Advanced Satellite for Cosmology and Astrophysics 
(Tanaka, Inoue, & Holt 1994). 
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Soft X-rays can also arise from an elliptical's general interstellar medium (Fabian & Canizares 
1988; Mahadevan 1997; Di Matteo et al. 2000). For the galaxies in this study, a Bondi analysis 
of that medium can produce an estimate for the Bondi accretion rate, Mb, which can then be 
compared with the absolute ADAF accretion rate, Ma, derived at 8.5 GHz. The black hole masses 
from Magorrian et al. (1998), in combination with = 1 being typical for elliptical galaxies 
(Di Matteo et al. 2000), results in the Bondi radii listed in the table. Assuming further that the 
pressure, P = W^Pq cm~^ K, at the Bondi radius satisfies Pe = 1 — W (Di Matteo et al. 2000), then 
the table gives the corresponding range in Bondi accretion rates, Mb. The Bondi rate estimates 
for Pq = I are generally consistent with the limits on the absolute ADAF rates, Ma) imposed 
at 8.5 GHz, although there is an order-of-magnitude discrepancy in the case of NGC4291. In 
constrast, for Pq = 10 the Bondi rates exceed the ADAF rates for all four galaxies. Meaningful 
comparions between Mb and Ma must clearly await a reliable pressure profile for ellipticals in 
general or, better still, for the galaxies in this study. Futhermore, integration over any adopted 
pressure profile should not violate Lrass listed in the table. 

Further radio observations to quantify any jet contamination in NGC 4621 and NGC 4660 are 
also needed, for two reasons. First, removal of jet contamination leads to reduced ADAF emission, 
which in turn implies even lower ADAF accretion rates. Second, the presence of jet emission could 
signify the importance of physical processes ignored in canonical ADAF models and explored by 
Di Matteo et al. (2000). Follow-up VLA imaging at other frequencies, v, but at matched angular 
resolutions would help assess jet contamination: the flux density from an ADAF-dominated source 

1 2 

is expected to rise as i'3~5 (Mahadevan 1997; Yi &; Boughn 1998), whereas the flux density from 
a jet-dominated source is expected to exhibit a flatter spectral slope. Although the 500-;uJy ba~^ 
upper limits at 5 GHz from Wrobcl &; Heeschen (1991) do have similar resolutions, they are too 
insensitive to impose useful constraints on spectral slopes between 5 and 8.5 GHz. Moreover, higher- 
resolution imaging with the NRAO VLBA could provide morphological evidence for nonthermal 
synchrotron jets and spatially separate jet from ADAF emission, as was successfully done for the 
spiral galaxy NGC 4258 (Herrnstein et al. 1998). Finally, the image rms values, cts.sghz, achieved 
in this study required 2 hours of integration time. Similar values for cJg.sGHz will be achieved after 
a 2- minute integration with the expanded VLA (NRAO 2000). This will make it feasible, for the 
first time, to target all radio-quiescent ellipticals studied by Magorrian et al. (1998), which should 
significantly advance our understanding of advection-dominated accretion flows as a new class of 
extragalactic radio emitters. 

The authors thank Dr. R. Mahadevan for discussions. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, 
Caltcch, under contract with the National Aeronautics and Space Administration. The Digitized 
Sky Surveys were produced at the Space Telescope Science Institute under U.S. Government grant 
NAG W-2166. The National Geographic Society - Palomar Observatory Sky Atlas (POSS-I) was 
made by the California Institute of Technology with grants from the National Geographic Society. 
The Second Palomar Observatory Sky Survey (POSS-II) was made by the California Institute of 
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Table 1. Quiescent Elliptical Galaxies at 8.5 GHz 



Symbol 


Units 


NGC4291 


NGC 4564 


NGC 4621 


NGC 4660 


Reference 






2.0 


2.6 


4.5 


2.4 


1 


fDSS 


// 


1.7 


2.0 


3.0 


1.8 


2 


-Rs.SGHz 


// // o 


15.4,8.4,+56 


10.5,9.0,-22 


10.2,9.0,4-19 


14.7,9.1,-^54 


3 


C^S.SGHz 


/iJy ba~^ 


16 


13 


13 


14 


3 




/xJy ba~^ 


< 64 


< 52 


153±14 


143±15 


3 


D 


Mpc 


28.6 


15.3 


15.3 


15.3 


4 


s 


PC/" 


139 


74 


74 


74 


4 


^8.5GHz 


W Hz-i 


< 6.3 X 10^8 


< 1.5 X 10^8 


< 4.3 X lO^s 


< 4.0 X lO^s 


3 


ms 


lO^M© 


19 


2.5 


2.8 


2.8 


4 


rriA 




< 1.6 X 10-^ 


< 3.6 X 10-^ 


< 7.8 X 10-^ 


< 7.4 X 10-^ 


3 


Me 


M© yr-i 


42. 


5.5 


6.2 


6.2 


4 


Ma 


Mq yr-i 


< 0.007 


< 0.002 


< 0.005 


< 0.005 


3 


-^^RASS 


erg s^^ 


8.0 X 10"^° 


< 4.2 X 10^^ 


< 4.6 X 10^^ 


< 3.9 X 10^^ 


5 


-^A,lkeV 


erg s^^ 


< 5.9 X 10^8 


< 4.0 X 10^8 


< 2.1 X 10^9 


< 1.9 X 10^9 


3 


-^A,6keV 


erg s"^ 


< 3.5 X 10^9 


< 2.3 X 10^9 


< 1.2 X 10^° 


< 1.1 X 10^° 


3 


r-B 


pc 


82 


11 


12 


12 


4,6 


Mb 


Mq yr-i 


0.07-0.7 


0.001 - 0.01 


0.001 - 0.01 


0.001 - 0.01 


4,6 



References. — (1) Wrobel 1991; (2) Cotton et al. 1999; (3) This work; (4) Magorrian et al. 
1998; (5) Beuing et al. 1999; (6) Di Matteo et al. 2000. 
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Fig. 1. — VLA images of Stokes / emission from the inner ±1' of four elliptical galaxies at a 
frequency of 8.5 GHz. Contour levels are at multiples of the bottom contour level, which is twice 
the image rms value, crg.sGHz- Hatched ellipses show the restoring beam areas at FWHM. Negative 
contours are dashed and positive ones are solid. Upper left: NGC4291. Upper right: NGC4564. 
Lower left: NGC4621. Lower right: NGC4660. 
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